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I.  INTRODUCTION 


In  1958  Karpov  and  Bradley1  observed  short  flights  of  shell,  accompanied 
by  substantial  reductions  In  spin  rates.  An  eleven  percent  range  reduction, 
for  example,  was  accompanied  by  a  fifty  percent  reduction  In  spin  over  thirty 
seconds.  Later  tests  showed  that  the  spin  loss  occurs  after  the  coning  motion 
exceeds  ten  degrees  and  Is  not  the  cause  of  growth  In  the  coning  motion.  It 
was  shown2-3  in  1977  that  this  coning  growth  and  spin  decay  could  be  explained 
by  the  motion  of  loose  internal  components.  It  should  be  emphasized  that  In 
all  monitored  flights  of  projectiles  with  large  coning  motion  and  no  loose 
Internal  parts,  no  large  despin  moments  have  been  observed.  Thus,  the 
combined  presence  of  a  large  coning  motion  and  a  large  despin  moment  Is 
symptomatic  of  a  loose  Internal  component. 

In  the  summer  of  1974,  a  155mm  shell  with  a  liquid  payload  was  fired  with 
a  yawsonde  telemetry  unit4  and  this  projectile  developed  a  coning  motion  with 
amplitude  In  excess  of  forty  degrees.5  When  the  coning  motion  exceeded  forty 
degrees,  the  telemetry  record  showed  a  very  rapid  despin  of  forty  percent  In 
less  than  five  seconds.  Since  1974,  six  more  observations  of  liquid-filled 
155mm  projectiles  with  large  coning  motions  and  large  despin  rates  have  been 
made.6-8  In  addition,  thirty-six  observations  have  been  made  of  large  coning 


1.  B.G.  Karpov  and  J.V.  Bradley ,  "A  Study  of  Causes  of  Short  Ranges  of  the 
8 "  TS17  Shell ,"  BRL  Report  1049 ,  May  19S8  (AD  377S48 ). 

2.  C.H.  Murphy,  " Influence  of  Moving  Internal  Parts  on  Angular  Motion  of 
Spinning  Projectiles,"  Journal  of  Guidance  and  Control ,  Vol.  1,  March - 
April  1978,  pp.  117-122.  "See  also  bRL-MR-2731,  February  1977 

(AD  AO 3? 538). 

3.  V.G.  Soper ,  "Projectile  Instability  Produced  by  Internal  Friction,"  AIAA 
Journal.  Vol .  18,  Jan  1978,  pp.  8-11. 

4*  V.H.  Msrmgen,  "Measurements  of  the  Dynamical  Behavior  of  Projeotiles 
Over  Long  Flight  Paths,"  Journal  of  Spacecraft  and  Rockets,  Vol .  8,  April 
1971,  pp.  380-385.  See  also  BRL  MR-2079,  November  1970  (AD  717002). 

5.  V.P.  D'Amico,  V .  Oskay,  V.H.  Clay,  "Flight  Tests  of  the  155m  XM687 
Binary  Projectile  and  Associated  Design  Modification  Prior  to  the  Nioolet 
Vinter  Test  1974-1975,"  BRL-MR-2748,  May  1977  (AD  B019969). 

6.  V.P,  D'Amico,  V.H.  Clay,  and  A.  Mark,  "Yawsonde  Data  for  M687-Type 
Projeotiles  with  Application  to  Rapid  Spin  Decoy  and  Stewarteon-Type 
Spin-Up  Instabilities,"  ARBRL-MR-03027 ,  June  1980  (AD  A089646). 

7.  V.P.  D'Amico,  V.H.  Clay,  and  A.  Mark,  "Diagnostic  Tests  for  Vick-Type 
Payloads  and  High-Viscosity  liquids ,*  ARBRL-MR-02913 ,  April  1979 

(AD  A072812). 

8.  V.P.  D'Amico  and  V.  H.  clay,  "High  Viscosity  Liquid  Payload  Yawsonde  Data 
for  Small  Launch  Yaws,"  ARBRL-MR-03029,  June  1980  (AD  A088411). 
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notions  and  despin  moments  for  projectiles  with  payloads  of  liquid-solid 
mixtures;  namely,  twenty-six  XM61  projectiles  containing  liquid  plus  cotton 
wicks7*9*11  and  ten  XH825  projectiles  containing  liquid  plus  felt 
wedges.12*14  Indeed,  all  observed  flights  of  projectiles  with  liquid  or 
liquid-solid  payloads  that  performed  large  coning  motion  exhibited  large 
losses  In  spin.  .  V 

This  characteristic  association  of  a  large  despin  moment  with  a  large 
coning  motion  for  a  projectile  with  a  moving  payload  can  become  a  diagnostic 
tool.  Hiller15  suggested  using  this  tool  to  determine  the  existence  of  small 
amplitude  unstable  11 quid- Induced  side  moments.  Spin  measurements  made  during 
coning  motion  on  a  gyroscope  predicted  flight  pitch  instabilities  caused  by 
very  viscous  liquids,  and  these  were  observed  In  flight.16  The  original  M687 
flight  Instability  was  caused  by  a  large  external  Magnus  moment,  even  though 
the  liquid  caused  a  rapid  despin.  Thus  the  assumed  connection  between  large 
coning  motion  roll  moments  and  small  coning  motion  side  moments  may  have  been 
fortuitous  for  the  case  of  very  viscous  payloads. 

The  linear  liquid-induced  side  moment  was  first  computed  by  Stewartson17 
for  an  Invlscld  liquid  payload  by  use  of  elgenfrequencles  determined  by  the 


9,  V.P.  D^Amioo,  "Early  Flight  Experiences  Vith  the  XH781 ,"  BRL-MR-2791, 
September  1977  (AD  B024975L). 

10 ,  V.P.  D'Amico,  " Field  Taste  of  the  XH761:  First  Diagnostic  Test ,"  BRL-MR- 
2792 ,  September  1977  (AD  B024976L), 

11*  V.P.  D'Amico,  "Field  Teste  of  the  XM761:  Second  Diagnostic  Test,"  ABBRL- 
HR-02806 ,  January  1978  (AD  B025805L) , 

22,  V,P,  D'Amico ,  "Aeroballietio  Testing  of  the  XM825  Projeotile:  Phase  I," 
ARBRL-MR-02911,  March  1979  (AD  B0S7680). 

18,  V.P.  D'Amico,  "Aeroballietio  Testing  of  the  XN825  Projeotile:  Phase  II," 
ARBRL-MR-08072,  January  1981  (AD  AD98038), 

14,  V.P.  D'Amico,  "Aeroballietio  Testing  of  the  XH825  Projeotile:  Phase  III, 
High  Hues  Is  Velocity  and  High  Quadrant  Elevation,"  ARBRL-HR-03196 , 
September  1982  (AD  B088511L), 

18,  H,C,  Hiller,  "Flight  Instabilities  of  Spinning  Projectiles  Having  Hon - 
rigid  Payloads,"  Journal  of  Guidance,  Control,  and  Dynamics,  Vol,  8, 
March- April  1982,  pp,  251-187, 

16,  V.P.  D'Amico  and  H,C,  Hiller,  "Flight  Instability  Produoed  by  a  Rapidly 
Spinning,  Highly  Viscous  Liquid,"  Journal  of  Spaoeoraft  and  Rockets,  Vol, 
14,  Jan-Feb  1979,  pp,  62-84, 

17,  K,  Stevartson,  "On  the  Stability  of  a  Spinning  Top  Containing  liquid," 
Journal  of  Fluid  Mechanics,  Vol,  5,  Part  4,  September  1959,  pp,  577-592, 
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fineness  ratio  of  the  cylindrical  container.  Wedemeyer18  Introduced  boundary 
layers  on  the  walls  of  the  container  and  was  able  to  determine  viscous  correc¬ 
tions  for  Stewartson's  elgenfrequencles,  which  could  then  be  used  In  Stewart- 
son's  side  moment  calculation.  Murphy19  then  completed  the  linear  boundary 
layer  theory  by  Including  all  pressure  and  wall  shear  contributions  to  the 
liquid-induced  side  moment.  Tne  Stewartson-Wedemeyer  eigenvalue  calculations 
have  been  Improved  at  low  Reynolds  number  by  Kitchens  et  al20  through  the 
replacement  of  the  cylindrical  wall  boundary  approximation  by  a  linearized 
Navler-Stokes  approach.  Next, Gerber  et  al2^22  extended  this  linearized  NS 
technique  to  compute  better  side  moment  coefficients  for  Reynolds  numbers  less 
than  10,000. 

The  only  theoretical  work  on  liquid-induced  roll  moments  was  done  by 
Vaughn23  In  1978.  Although  fair  agreement  with  Miller's  data  was  obtained, 
the  work  was  marred  by  sane  hard-to-justlfy  algebraic  steps.  It  Is  the 
purpose  of  this  report  to  compute  the  roll  moment  associated  with  the  linear 
perturbation  flow  of  Reference  19. 


II.  LIQUID  ROLL  1  OMENT 

The  pitching  and  yawing  motion  of  a  symmetric  spinning  projectile  can  be 
represented  as  the  sum  of  two  rotating  exponentially  growing  two-dimensional 
vectors.  In  terms  of  complex  variables,  this  relationship  assumes  the  form214 


18,  B,N,  Wedemeyer,  "Viscous  Correction  to  Stewartaon'  a  Stability  Criterion." 
BRL-R-132S,  June  1966  (AD  489687). 

19,  C,H,  Murphy ,  "Angular  Motion  of  a  Spinning  Projectile  with  a  Viaoouo 
Liquid  Payload ,"  Memorandum  Report  ARBRL-MR-03194 ,  August  1982  (AD 
A116676).  See  also  Journal  of  Guidance,  Control,  and  Dynamics,  Vol,  6, 
July- August  1983 ,  pp,  280-286, 

20,  C.U,  Kitchens ,  Jr.,  N,  Gerber ,  and  R,  Sedney,  "Oscillations  of  a  Liquid 
in  a  Rotating  Cylinder:  Solid  Body  Rotation,"  URL  Technical  Report  ARBRL- 
TR-02081 ,  June  1978  (AD  A057759). 

21,  8.  Gerber,  R,  Sedney,  and  J.M,  Sartos,  "Pressure  Moment  on  a  Liquid- 
Filled  Projsotile:  Solid  Body  Rotation,"  BRL  Technical  Report  ARBRL-TR - 
02422,  Ootober  1982  (AD  A12QS67 ). 

22,  N.  Gerber  and  R,  Sedney,  "Moment  on  a  Liquid-Filled  Spinning  and  Nutating 
Projsotile:  Solid  Body  Rotation,"  BRL  Technical  Report  ARBRL-TR-02470 , 
February  1983  (AD  A125332). 

23,  H,R,  Vaughn,  " Flight  Dynamic  Instabilities  of  Fluid-Filled  Projectiles ," 
Sandia  Laboratories,  Albuquerque,  NM,  SAND  78-0999,  June  1978. 

24,  C,H,  Murphy,  "Free  Flight  Motion  of  Symmetric  Missiles,"  BRL-R-1216, 

July  1983  (AD  442757). 


(2.1) 


c  si*  1  a  -  Kje1*1  *  ty*2 


where  In  (Kj/Kj0)  •  Tj|+|t 


♦j  *  ♦jO  *  *nd 

+  Is  the  roll  rate. 

The  transverse  moment  exerted  by  a  liquid  payload  was  expressed  In 
Reference  19*  In  the  following  form. 


where 


♦  1  \t  "  \  »  ♦  t'l  CLM1  V  *  ^  *2e  3  (2. 


*  mass  of  liquid  In  fully  filled  container 


a  *  maximum  radius  of  liquid  container. 


For  linear  fluid  motion,  should  depend  on  tj,  cj,  time,  Reynolds  number, 

fill  ratio,  shape  of  cavity,  and  direction  of  spin.  A  similar  remark  applies 
t0 

The  C«N  are  complex  quantities  whose  Imaginary  parts  represent  In-plane 

j 

moments  causing  rotation  In  the  plane  of  exp  (1+j)  and  whose  real  parts 
represent  side  moments  causing  rotation  out  of  the  plane  of  exp  (1^).  Thus, 
Cjj|  with  Its  dependence  on  the  direction  of  spin  can  be  written  as 

J 


:LMj  "  Y  CtSMj  *  1  ClIMj 


(2.3) 


where  ,  CLJH  are  real  and  y  •  i/|  j|. 

J  j 


*Sinae  relations  from  Reference  19  ore  used  throughout  this  report,  on  errata 
for  Reference  19  is  given  in  Appendix  4. 
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The  liquid  roll  moment  can  now  be  defined  In  a  similar  way.  Symmetry 
considerations  Imply  an  even  dependence  on  coning  amplitudes  and  an  odd 
dependence  on  spin  rate. 


*’•  "  ”1  a  ^LRMq  *  T1  K1  CLRM,  +  t2  K2C  CLRM^ 


’1 


(2.4) 


where  C^u  are  functions  of  Reynolds  number,  fill  ratio,  shape  of  cavity,  and 

time.  In  addition,  for  j  ■  1  and  2,  the  may  be  functions  of  Tj,  and 

Kj  •  is  roll  moment  coefficient  exerted  by  the  liquid  during  spin- 

up  and  Is  zero  when  steady  state  Is  reached.  For  the  steady-state  coning 
motion  considered  In  the  Stewartson-Uedemeyer  theory  (K*  ■  0),  only  one  roll 

moment  coefficient  Is  present  and  we  will  omit  the  subscript  on  tj  and  : 


t  Rj  ^LRM* 


(2.5) 


Throughout  this  report,  we  will  consider  only  roll  moments  that  can  be 
approximated  by  the  one-term  Eq.  (2.5)  rather  than  thq  more  complete  three- 
term  Eq.  (2.4)  and  will  consider  only  positive  spin,  *. 

The  roll  rate  of  a  projectile  In  flight  Is  controlled  by  the  sum  of  the 
internal  liquid  roll  moment  and  the  external  aerodynamic  roll  damping 
moment.  For  pure  coning  motion 

lx  ♦  ■  (-^-)  St  (^p)  ♦  m^  a^  &  t  (2.6) 

For  the  liquid  payloads,  the  large  Increase  in  roll  moment  usually  occurs 
suddenly  when  the  coning  motion  amplitude  exceeds  30*.  On  the  other  hand,  the 
liquid/solid  XM761/XM825  payloads  produce  large  roll  moments  at  a  much  lower 
coning  amplitude  of  10-15*.  For  all  payloads,  the  coning  motion  Increases 
rapidly  to  a  steady-state  amplitude  of  40-45*. 

Since  definitions  (2.1,  2.2,  2.4)  were  for  small  coning  angles,  more 
precise  expressions  are  required  for  angles  In  excess  of  ten  degrees.  The  yaw 
plane  is  the  plane  containing  the  missile  axis  and  the  velocity  vector.  The 
angle  between  the  two  vectors  is  the  total  angle  of  attack,  ay.  The  amplitude 

of  the  complex  variable  |  is  sin  <y  and  the  argument  of  |  is  the  phase  angle 

of  the  yew  plane  about  the  missile  axis.  For  the  pure  coning  cation  of  Eq* 
(2.S),  the  exact  definition  of  Kj  is 


■  III  ■  sin  ^ 


(2.7) 
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Unfortunately,  the  yawsonde  data  for  the  XM761/XM825  flights  were  too 
poor  to  yield  numerical  values  of  the  liquid  roll  moment  coefficient.  Seven 
flights  of  the  N687  payloads  did  have  sufficiently  good  data  and  their  liquid 

roll  moment  values  could  be  computed  from  Eq.  (2.6).  The  change  In  roll 

•• 

acceleration  at  high  cone  angles,  A  4,  can  be  easily  obtained  from  the  yaw* 
sonde  data.  From  Eq.  (2.6) 


4  ♦  -  <\  *X>  CLRH 


(2.8) 


Table  1  gives  the  seven  M687  values  of  the  roll  moment  coefficient.  For 
very  low  Reynolds  number,  a  value  of  about  -.04  seems  representative  while  a 
high  Reynolds  number  value  of  -.01  Is  Indicated.  Since  Round  7254  was  not 
completely  filled.  It  Is  an  Isolated  value.  Miller's  gyroscope  data15  1$ 
given  In  Figure  1  and  Indicates  a  range  of  -.001  to  -.05  for  CLRM.  Two 

recently  obtained  values  of  roll  moment  for  eight-inch  shell25  are  also 

31ven  In  Table  1.  Since  the  flight  data  were  poorly  resolved  during  spin- 
own,  the  liquid  roll  moment  coefficients  of  Table  1  should  be  considered 
estimates  with  accuracies  no  better  than  30X. 


III.  LIQUID  ROLL  MOMENT  THEORY 

Two  coordinate  systems  will  be  used  In  this  report:  the  nonrolling 

aeroballlstlc  XYZ  system  whose  X-axis  1$  fixed  along  the  missile's  axis  of 
symmetry  and  the  Inertial  XYZ  system  whose  X-axis  Is  tangent  to  the  trajectory 
at  time  zero.*  Both  coordinate  systems  have  origins  at  the  center  of  the 
cylindrical  payload  cavity,  which  Is  assumed  to  be  at  the  center  of  mass  of 
the  projectile.  Location  In  the  cavity  can  be  specified  in  the  aeroballlstlc 

system  hy  the  cylindrical  coordinates  x,  r,  e  and  In  the  Inertial  system  by  x, 

r,  a.  The  boundary  of  the  cavity  Is  given  by  x  •  $c  and  r  ■  a  where  2c  is  the 
height  of  the  cavity.  Linear  relations  between  cylindrical  coordinates  In  the 
two  coordinate  system*  were  derived  In  Reference  19: 

x  •  x  -  r  Kj  cos  -  •)  (3.1) 


r  »  r  ♦  x  cos  (4j  -  a)  (3.2) 


*  The  trajectory  sok  he  well  approximate  by  a  straight  line  for  a  • umber  of 
period*  of  the  angular  motion, 

26.  IT.  P,  O' Mtioo  and  8.  «/.  talamnohili,  “taooonlo  feete  of  the  JtVJ?? 

Binary  Pro j eat  He:  Photo  BRL  Memorandum  Popart  in  preparation . 
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TABU  l.  EXPERIMENTAL  VALUES  OF  C 


Fill  r#Uo  wit  .87  for  Round  7254;  all  other  rounds  were  1001  fill 
*•  8- Inch  XM877  projectiles. 


-  •  \»/’r)  Kj  sin  Uj  -  8)  j3  3j 

the  fn'thVSir^tlln'ofMcfeesInS 


|~  .  /£.£  *!.  *  /3f\2  ^  /3f  \21'1/Z  , 

X  \8x*  ap*  r36/^V3x/  +  It rj  +  (rje/  J  =  (e«>  exr*  ex8>  (3.4) 


*x  s  <exx»  exr»  exe> 


*  (l,  -  ^  cos  -  e),  -  Kj  sin  -  e)  j 


(3.5) 


%  =  ^erx*  err*  er0) 


=  cos  -  a),  i,  Kj  (x/r)  sin  -  e)\ 


(3.6) 


%  =  (e3x*  e9r*  ®?e) 

*  (Ki  sin  (♦j  -  8),  -  Kj  (x/r)  sin  -  0),  1) 


(3.7) 


Eqs.  (3.5  -  3.7)  are  correct  to  terms  of  order  Kj2. 

..  S^nce  c°mPr®ss1ve  viscous  shear  Is  neglected  in  boundary  layer  flows 
veryCSp?rfom26the  shear  tensor  in  cy11ndr'cal  coordinates  have  the 


Txx  *  Trr  *  Tee  *  -P 


T  ■  T 

xr  Vx 


J"'.  !!*! 

L**  3r  J 


(3.8) 


(3.9) 


48  ‘  nM"  XoOicu-Hill  Book  Company,  Bow 


KEVzKSSSSXSl 


^ w"^Tvv-T  V  ^Tvv?vv! 


v 


Txe  3  Tax 


(3.10) 


(3.11) 


The  stress  vector  on  an  area  element  of 
surface  of  the  liquid  container  is 


the  Interior  lateral 


cylindrical 


dF*9> 


where 


(3.12) 


dFtx  3  •  [efx  *xx  +  err  Trx  +  «?8  *3x1  *  d*  d» 

J  r=a 

dFtr  =  *  [erx  Txr  +  err  Vr  *  e?9  V  1  a  d* 

J  cw 

r«a 

dF*e  ’  -  [erx  *xe  +  err  're  +  ere  'eej.  a  d*~  d« 

r*a 

The  corresponding  roll  moment  differential  is 
dMxi  a  a  h  ’  d(^ 

,  /  (3.13) 

‘  ”  l  Tre  +  CV  sin  ^i  "  ®)  +  cos  (^  -  9)]  kx  l  a2  dx  dS 

r«a 

SllLflKS.JSCy"  ‘he  'ater41  d"»  normal  to  the 


SV> 

Txr  "  w  aF 


TX0*  0 


(3.14) 

(3.15) 
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Tre  =  u  r 


»(yr) 

5F 


(3.16) 


The  roll  moment  differential  on  the  forward  and  rear  end  walls  (j  =1,  2) 
can  be  computed  in  a  similar  manner. 


<xej  3  M^e-  Cv  <*/p)  s1n  ^r0) 


(3.17) 


+  Tr0  cos 


(<t*re)]  U 


,A»  ,/V 

r  dr  de 


On  the  endwall  only  derivatives  with  respect  to  x  are  important. 

3  V 

r 

Txr  v  3x 


Txe  u  3x 


Tre  3  ® 


(3.18) 

(3.19) 

(3.20) 


We  will  now  compute  the  quadratic  liquid  roll  moment  under  the  restric¬ 
tive  assumption  that  only  linear  fluid  mechanics  terms  need  be  considered.  In 
Reference  19,  the  velocity  components  were  expanded  in  terms  of  perturbation 

A 

velocities  that  were  linear  in  K. 


-  1  8  |  (,}) 

(3.21) 

(3.22) 

{"s  es*  '  1e}  (aj) 

(3.23) 

where 


s  ■  (e  +  i)  t 
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trjXt 


**  \.  x  -  V  -  *  .  *.« 


and  where  u  ,  v  ,  w  are  functions  of  r  and  x  and  are  linear  In  k.  Eqs. 

3  J  5 

(3.21-3.23)  give  the  velocity  components  In  terms  of  x,  r,  6.  In  computing 
the  linear  liquid  side  and  In-plane  moments,  x,  r,  0  could  simply  be  replaced 

by  x,  r,  e.  The  liquid  roll  moment,  however.  Is  a  quadratic  function  of 

and  Eqs.  (3. 1-3.3)  must  be  used.  A  useful  approximation  for  this  change  of 
variables  is 


h  (r»x)  e-ie  =  h  (r,x)  e"1^  +  (r  -  r) 


+  (x  -  x)  -  1h  (r,x)  (0-0)  e'ie 

The  roll  moment  induced  by  the  lateral  wall  can  now  be  computed  by 
Integrating  Eq.  (3.13)  with  the  aid  of  Eqs.  (3.14-3.16,  3.21-3.24). 


(3.24) 


r  f  m  *2*2  v  2 \— 1 

dm  =  ( t  m,  a  (j>  K.  ) 


c  >2ir 


u 


=  (2tC  )• 


n 


(3.25) 


where 


f  =  ax 
£ 


K  •  K10  e 


.2 


Similarly,  Eq.  (3.17)  can  be  integrated  to  yield 
CLRMel  ‘  CLRMe2  s  <1/2>  CLRMe 


■  (2Tac)’5  (  R  i  f e  K’1!  r  dr 
JO  I  )r-r 


(3.26) 


(3.27) 


(3.28) 


2  As  a  /  3  ws 

where  f0  -  -r4  — ^  +  c  4s-  [  w.  +  1  vt  r  r - 

e  ax4  ax  i  s  s 


■  k  ■  •  *  \  »  ■ 


,  /.  4\  4\  %  A  %  %  A  A  A 


(3.29) 


cnitsxs: 


Finally 


CLRM  a  CLRM.  +  CLRM# 
e  t 


In  Appendix  B  it  is  shown  that 

CLRM  a  "CLSM  +  (Te/2)tl 


(4/3)  (c/a)2] 


(3.30) 


(3.31) 


This  equation  is  valid  for  linear  fluid  mechanics  and  states  that  for  the 
simple  case  of  constant  amplitude  motion  (e  *  0)  the  liquid  roll  moment  is  the 
negative  of  the  liquid  side  moment.  If  the  complete  nonlinear  fluid  equations 
are  used,  additional  contributions  to  the  quadratic  liquid  roll  moment  which 
are  large  enough  to  affect  Eq.  (3.31)  may  be  present. 


IV.  DISCUSSION 

Reference  2  considers  the  effect  of  two  types  of  internal  component 
motion:  (a)  a  mass  moving  in  a  circle  In  a  plane  normal  to  the  missile's 
axis;  (b)  forced  coning  motion  of  a  spinning  component.  For  these  motions  the 
side  moment,  M$,  and  tne  roll  moment,  M^,  were  computed  and  for  both  cases 


"sr*  -Ki 


(4.1) 


This  Is  precisely  Eq.  (3.31)  for  constant  amplitude  motion.  Equation  (4.1)  is 
equivalent  to  the  statement  that  the  liquid  moment  about  the  trajectory  Is 
zero,  l.e., 

Mu  «  0  .  (4.2) 

This  essentially  follows  from  the  linear  assumption  that  the  azimuthal 
velocity  about  the  trajectory  varies  as  exp  (-18).  A  nonlinear  theory  could 
well  have  terms  that  are  independent  of  8. 

The  liquid  side  moment  was  computed  in  Reference  19  for  c/a  =»  4.291  and 
Re  from  102  to  106;  the  results  are  given  In  Figures  2-6.  Miller15  observed 
the  roll  moment  to  be  Independent  of  roll  rate.  This  means  that  the  liquid 
roll  moment  should  be  proportional  to  t.  As  can  be  seen  from  Figures  2-6, 
this  Is  roughly  the  case  for  the  liquid  side  moment  coefficient.  Since 
Miller's  measurements  were  taken  for  t  varying  from  .12  to  .25,  we  will  use 
the  liquid  side  moment  coefficient  for  t  *  .18  to  compare  with  Miller's 
measurements  In  Figure  1.  The  agreement  Is  encouraging  when  we  recall  that  we 
are  applying  a  linear  fluid  mechanics  calculation  to  measurements  taken  at  a 
20°  cone  angle  and  a  boundary  theory  to  Reynolds  numbers  as  low  as  100. 
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V.  SUMMARY 

A  .  \\  Ih®  !1neaC  ttieor>  for  a  coning  and  spinning  liquid  payload  has  been 
extended  to  the  calculation  of  a  quadratic  roll  moment. 

exact thier.rn!fJJlnLq!kdr?Jtc  ro,!Jmoment  coefficient  is  shown  to  be 
?9at1ve  °f  the  1inear  side  moment  coefficient  for  constant 
amplitude  coning  motion. 

nhc  Jh1s  relationship  between  side  moment  and  roll  moment  has  been 
supped  bj  Pe;MtSfSM,™ld"?vfJlds?0n’p0"entS  and  is  reasonab^ 


DESPIN  DATA  BY  MILLER 


4 
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APPENDIX  A 

ERRATA  FOR  REFERENCE  19 

Third  term  of  Eq.  (5.20)  should  contain  (s  -  1),  not  (s  - 
Second  line  of  Eq.  (6.4)  should  read: 


1  (c/2a) 


j  ;  [C?]r-a  *M  *  +  (h/C)?  "pth 


3.  Eq.  (7.1)  should  read: 


%  ■  1  *  "e)'1  e's*  f  ‘  f2V8  «;t  de  dx' 

J-l  Jo 


4.  Eq.  (7.4)  should  read: 


(6.4) 


(7.1) 


**e  '  (a  *  Re)  1  (  f  *  It  -  '•s.)  1 .  r<lr  *  (h/c)2  B  (7.4) 

L  Jx  •  -1 


where 


»u.h  "  <2  *  **>’1  U/»h) 


h  «■  j;.  A 


5.  Third  term  of  Eq,  (9.4)  should  contain  (s  -  i)t  not  (s-1), 
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APPENDIX  B 

Relationship  Between  CLRM  and  CL$M 


APPENDIX  B 

RELATIONSHIP  BETWEEN  CLRM  AND  CL$M 


In  Reference  19 
lateral  and  endwalls 
trlbutlons  are  given 
complete  side  moment 


SrefSU4e  Parts.°f  the  side  moment  coefficient  for  the 

ineEasVe?7  5  E7Sa^6*tk’  6*5^**  The  visc<>us  shear  con- 
jlILSfV7'2*  7*4]«  These  can  be  combined  to  yield  the 
coefficients  for  the  lateral  and  endwalls,  respectively: 


(B-l) 


(B-2) 


(B-3) 


■  *  1  Ps  +  c  Tx  (ws  -  1  *s)  Re'1 


(B-4) 


As  Is  shown  by  Eqs. 
very  similar  form: 


(3,25  -  3.29),  the  roll  moment  coefficients  have  a 


(B-5) 


x«c 


(B-6) 


*  Brackets  identify  equations  from  Reference  19, 
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R3W5 

I 


I 


1 


I 

R 


3 

.’5 

■X 


where 


V 


2 

3  W. 


Since 


L 

3r2 

-r2 

3S 

- 

"17 

2 

3ws 

(a»x) 

2 

-a2  — i  .  i,_i 


3r3x 


3ll( 

3r 


Re 


-1 


(B-7) 


+  c  r 


3  w$  aw 
3xar  ax 


s  +H+  1  ”* 


3x 


Re' 


(B-8) 


3r 


2  is  the  ^Inant  term  In  CLRM  ,  It  must  be  computed  more 

l 


ferentlal  ^  U?Jal  boun ^  uy*r  approximation.  The  linear  circum 

rerentlal  momentum  equation  as  given  by  Eq.  [B6J  is 


a2  Re 


Llr2  r3r  17~  "7  "  r2  J 


lap. 


(B-9) 


(s  -  1)  ws  +  2v$ - J. 


On  the  boundaries  r  =  a.  x  =  c  or  v  -  r  „  ***  .. 
conditions  as  given  by  Equations  [3.8  -  3’.9]S:  S  "“St  the  b°U"dar* 


WS  s  1  (s-1)  (x/a)  K  =  -1v 


(B-10) 


Thus,  on  the  boundary  surfaces  of  the  liquid  contal 


ner 


i.*s  >\  .2 

-  4.  -  +  - -  a  a 


3r 


r3r  ax 


r  iaps  2 

+  1  (s^  +  1)  (x/a)  K  |Re 


A 

;  f 


(B- 11) 


32w. 


On  the  lateral  boundary  (r  *  a) ,  — is  zero 


32W. 


b*  a"< 


2  13  teru  and  *~y  in  Eq.  (B-7)  can 


3r 


CLRM4  a  *  cLSMa  *  (2/3)  (c/a)2  xe 


(B- 12) 
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vL  *:-  *  ^*-*-  '  * 


2  2 

3  ws  ^S  -  3  ws 

Similarly,  — r-  and  — 2.  are  zero  on  the  endwalls  (x  =  ±c)  and  — a2,  in  Eq.  (B-8) 

3r  a r  3x 

can  be  eliminated  by  use  of  Eq.  (B-ll) 


+  te/2 


+  (t  Re) 


awR  . 

•-kk 


(B-13) 


The  derivative  of  w$  normal  to  the  endwall  at  the  junction  of  the  walls  can  be 

computed  from  Eq.  (8-10)  and  is  the  same  size  as  the  derivatives  omitted  in  the 
boundary  layer  approximations  of  Eqs.  (3.18-3.20).  Thus,  the  last  term  of  Eq. 
(B-13)  can  be  neglected.  Eqs.  (B-12  -  B-13)  can  then  be  added  to  yield  the  com¬ 
plete  liquid  roll  moment  coefficient: 


^LRM  55  ~^L$M  +  (Te/2)  [1  ~  (V3)  (c/a)  ] 


(B-14) 


For  constant  amplitude  coning  motion,  then,  the  liquid  roll  moment  is  the 
negative  of  the  liquid  side  moment. 
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APPENDIX  C 

LIQUID  ROLL  MOMENT  FOR  PARTIAL  FILL 

If  the  cylindrical  container  is  only  partially  filled  with  liquid,  a 
fully  spun-up  liquid  will  fill  the  space  between  the  outer  cylindrical  wall 
and  an  Inner  cylindrical  free  surface  with  radius  b.  The  ratio  of  the  volume 

2  2 

of  this  inner  cylinder  to  the  volume  of  the  complete  payload  cavity  is  b  /a  . 

2  2 

The  fill  ratio  for  the  payload  Is,  therefore,  1  -  b  /a  and  Is  denoted  by  f. 

In  the  presence  of  small  amplitude  coning  motion,  this  cylindrical  free 
surface  Is  deformed  slightly.  The  equation  for  this  surface  In  aeroballlstlc 


coordinates  Is 

r  =  b  (1  +  n)  (C-l) 

where  n  =  R  jns  (x)  es*  ‘  1e|  (C-2) 

A 

and  the  perturbation  function  n$  Is  linear  In  K. 

In  earth- fixed  variables  the  surface  has  the  form: 

F  (x,  r,  e,  t)  *  x2  +  b2  [1  +  n]2  -  x2  -  r2  *  0  (C-3) 

From  Eq.  [A12] 

x  »  x  -  R  |r  K  es*“  ie|,  (C-4) 


and  n  must  vary  such  that  the  free  surface  moves  with  the  liquid,  i.e.. 


DF 

ur 


v  + 

*r  3r  r  ae 


(C-5) 


Equations  (C-3)  -  (C-5)  with  Equations  (3.21  -  3.23)  can  now  be  used  to  obtain  an 

A 

equation  for  that  Is  linear  in  K. 


bhj  (x) 


-  x  K  + 


avs(b,x) 
s  -  1 


(C-6) 


According  to  Eq.  [B14],  at  x  «  *c 
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v$  ■  M*  -  0  (c/a)  K 


•  \  (c)  »  (-c)  «  0 


On  the  free  surface,  the  pressure  is  constant: 


(C-7) 

(C-8) 


*81 


(C-9) 


The  linear  version  of  this  Is  Eq.  [5.8],  which  was  used  as  the  Inner  boundary 
condition: 


(s  -  1)  ps  (b,x)  +  (b/a)  vs  (b,x)  *  0 


(C-10) 


The  side  moment  coefficient  has  the  same  Integral  relation  for  C,cu  as 
Eq.  (B-U  but  Eq.  (B-2)  becomes  l5"t 


:lsh .*  <t*c>1  0  *  {*•  *  "V?  ? 

•'b  x*c 


dr 


(C-ll) 


Similarly,  for  the  liquid  roll  moment  coefficient  Eq.  (B-5) Is  unaffected  but 
Eq.  (B-6) 


becomes 


C  -  («cr*  (*»{%**%? 

e  '»*  x-c 


Eq.  (B-13)  then  becomes 


*  *clsh  +  ( tc/2).  [1  -  (b/a)*] 


(C-U) 


(C43) 


•  ^LRH  *  "^LSH  *  Ci  •  (4/3)  (c/a)^  -  (b/a)*] 


(C— 14) 
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LIST  OF  SYMBOLS 


a 

c 


'LIN, 


'LM, 


radius  of  a  right-circular  cylindrical  cavity  containing  liquid 
liquid*19^  °f  *  r19ht*c1rcular  cylindrical  cavity  containing 

(j"2)  wode  11(Iu1d  1"-Plane  moment 
coefficients;  Imaginary  part  of  C,  M  . 

LMj 

and  slow  (j"2)  mode  complex  liquid  moment 
coefficients  defined  by  Eq.  (2.2) 

aerodynamic  roll  damping  moment  coefficient  defined  by  Eq.  (2.6) 
^LRMj  for  steady-state  coning  motion  (k2  •  0) 

Xid'Jy  E*  (M)(j"2)  ,,,,u1d  ro"  ““nt  eoerrie1«»« 

that  part  of  CLRM  Induced  by  the  end  walls  of  the  cavity 

that  part  of  CLRM  Induced  by  the  lateral  wall  of  the  cavity 

CLSMj  for  steady-state  coning  motion  (K2  •  0) 

real  part*ofC  BOde  liqu1d  s1de  ooaent  coefficients; 

***  *r»  ®1  un1t  vectors  along  the  aeroballlstic  cylindrical  axes 
*xx,exr,#xe  components  of  in  the  inertial  x,  r,  a  system 

erx,err**ra  c<*P«*nts  of  e?  In  the  Inertial  x,  r,  e  system 


'LRM 


'LRN, 


'LRM. 


'LRM. 


'LSM 


'LSM. 


e8**eV*e80  C08|PO«®nts  of  e^  In  the  Inertial  x,  r,  •  system 
*x  «1«1  moment  of  Inertia  of  the  projectile 

*  KI0®VO#i0) 


Kj  Magnitude  of  t)ie  fast  (j-1)  or  slow  (j»2)  yaw  arm 

Kjq  Initial  value  of  Kj 

t  reference  length 

•l_  mass  of  the  liquid  In  a  fully  filled  cylindrical  cavity 

H  »,  H  »,  H  »  XV2  components  of  the  moment  exerted  by  a  liquid  payload. 

LX  LY  12 

r,  r  radial  coordinates  In  the  Inertial  and  aeroballlstlc 

systems,  respectively 

Re  Reynolds  number 

s  (eel)  t 

S  reference  area 

t  time 

u.,  v_,  w_  components  of  the  liquid  velocity  perturbation  In  the 
555  Inertial  x,  r,  e  system 

V  magnitude  of  the  missile's  velocity 

vx.  Vr,  V#  velocity  components  In  the  Inertial  x,  r,  e  system 

x,  x  axial  coordinates  In  the  Inertial  and  aeroballlstlc  systems, 

respectively 

XV2  Inertial  exes;  X-axis  tangent  to  the  trajectory  at  time 

zero. 

XVZ  nonrolling  aeroballlstlc  axes;  Lexis  fixed  along  the 

missile's  axis  of  symmetry 

«•  1  angles  of  attack  and  side-slip  In  the  XYZ  system 

Y  I/I  II 

tj  nondlmenslonallzed  growth  rate  of  kj 

4,  I  azimuthal  coordinates  In  the  Inertial  and  aeroballlstlc 

systems,  respectively 

u  viscosity 
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€  ?♦  1  3 

p  air  density 

T  Tj  for  steady-state  coning  motion  (k2  ■  0) 

T<3  V**  the  noftdt«*ns1ona11zed  frequency  of  the  j-th  yaw  arm 

0-1.2) 

♦  roll  rate 

♦j  *J0  ♦  »j  *  ‘ 

♦jo  Initial  orientation  angle  of  the  j-th  yen  arm  (j-1,2) 
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